THE QUATERNIONIC DETERMINANT^ 
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Abstract. The determinant for complex matrices cannot be extended to quaternionic matrices. 
Instead, the Study determinant and the closely related (j-determinant are widely used. We show that 
the Study determinant can be characterized as the unique functional which extends the absolute value 
of the complex determinant and discuss its spectral and linear algebraic aspects. 
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1. Introduction. Quaternionic linear algebra is attracting growing interest in 
theoretical physics mainly in the context of quantum mechanics and field the- 

ory |(|. Quaternionic mathematical structures have recently appeared in discussing 
eigenvalue equations j?], |§| , group theory || [To) and grand unification model (ll], [I2] 
within a quaternionic formulation of quantum physics. 

The question of extending the determinant from complex to quaternionic ma- 
trices has been considered in the physical literature @]-@. The possibility of such 
an extension has been contemplated by Cayley E| , without much success, as early 
as 1845. A canonical determinant functional was introduced by Study jl4| and its 
properties axiomatized by Dieudonne Jlq| . The details can be found in the excellent 
survey paper of Aslasken Jl6| . Study's determinant is denoted as Sdet, and up to a 
trivial power factor, is identical to the g-determinant, det g , found in most of the lit- 
erature jlTj and to Dieudonne's determinant, denoted as Ddet. Study's determinant 
is closely related to the q-determinant and to Dieudonne's determinant. Specifically, 
det 9 = Sdet 2 = Ddet 4 . 

In these works, Sdet was considered as a generalization of the determinant, det, 
in the sense that the two functionals share a common set of axioms. Specifically, Sdet 
is the unique, up to a trivial power factor, functional T : M 71 x 71 which satisfies the 
following three axioms: 

1. J~{A) — if and only if A is singular; 

2. multiplicativity: T(AB) = T{A)F{B); 

3. T(l + rE i:j ) = 1 for i ^ j and reH; 

see JT6| . However, Sdet does not truly extend det. Indeed, the two functionals do not 
coincide on complex matrices, since the former is nonnegative while the latter is truly 
complex. In this paper we show that Sdet does extend the nonnegative functional 
det |, namely the two functionals coincide for complex matrices. More precisely, we 
show the following: 
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1) There exists no multiplicative functional on quaternionic matrices which coin- 
cides with dot on complex matrices. 

2) Sdet is the only non-constant multiplicative functional which coincides with 
det on complex matrices ( we remark that just like det[M] ^ 0, the inequality 
det[M]| ^ characterizes non-singular matrices over the complex numbers. The 
same central role in group theory over the quaternions will be played by Sdet[M] ^ 0). 

3) We show that the identities |dct(Af) = Y[ , m terms of eigenvalues, and 
det(M)| = Y[ °~i , in terms of singular values, extend to M quaternionic. Thus, al- 
though Sdet is originally defined through complexification (l^]-|2^], it can be given 
concrete spectral and numerical-analytic interpretations which do not require com- 
plexification. 

4) We show that the Schur complements identity for complex matrices, 

det [A] | |det[£> - CA~ l B]\ , 

extends to quaternionic matrices. 

5) We discuss formulas for det[ii], Adj[-ff] and based on the classical per- 
mutation and minor calculation (some of this material can be found in p2| ). It is 
interesting that this approach, pursued by Cayley without success in the context of 
general quaternionic matrices, is valid in the hermitian case. The functional det, de- 
fined this way for hermitian quaternionic matrices, is not multiplicative. Note that 
under the definition Sdet[Af] := det [M+M] one can extend the Study determinant to 
non-square matrices. 

In the last section we also discuss some open problems concerning the behavior of 
the determinant and the difficulties of extending the formula M -1 = Adj(M)/det(M) 
to quaternions. 

2. Notation. Quaternions, introduced by Hamilton [^3[ |24| in 1843, can be rep- 
resented by four real quantities 

q = a + ib + jc + kd, a, b, c, d G M , 

and three imaginary units i, j, k satisfying 

i 2 = j = k 2 = ijk = — 1 . 

We will denote by 

Re[g] := a and Im[<7] := q — a = ib + j c + k d , 

the real and imaginary parts of q. The quaternion skew-field H is an associative but 
non-commutative algebra of rank 4 over R, endowed with an involutory operation, 
called quaternionic conjugation, 

q = a — ib — jc— kd = Kc[q] — Im[q] , 

satisfying pq — qp for all q,p £ H. The quaternion norm \q\ is defined by 

\q\ 2 =qq = a 2 + b 2 + c 2 + d 2 . 
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Among the properties of the norm, to be used in subsequent sections, we mention 
here the following 



\pq\ = \qp\ = \q\ \p\ and |1 - pq\ = |1 - qp\ . 
Every nonzero quaternion q has a unique inverse 

q- 1 = q/\q\ 2 ■ 

Two quaternions p and q are called similar if 

q = s~ 1 ps , for some s£i. 

By replacing s by u = s/\s\, we may always assume s to be unitary. The usual 
complex conjugation in C may be obtained by choosing s = j or s = k. A necessary 
and sufficient condition for the similarity of p and q is given by 

Re [q] — Re[p] and Im[q]| = |Im[p]| . 

An equivalent condition is Re[q] — Re[p] and \q\ = \p\. Every similarity class contains 
a complex number, unique up to conjugation. Namely, every quaternion q is similar 
to Re[q] ± i \Im[q]\ . In particular, q and q are similar. It can be seen that s G H 
conjugates q and q (i.e. q — s~ 1 qs) if and only if Im[q] — or Re[gs] = Re[s] = 0. 
However, there exists no fixed s£l which conjugates q and q for all q £ H. 

3. Spectral theory. Spectral theory for complex matrices admits several possi- 
ble quaternionic extensions, which do not necessarily respect the fundamental theorem 
of algebra |?], ||, [2fJ-||29 . We shall be interested in the extension usually described 
as "right eigenvalues" j0T|], p9[ . 

Every n x n quaternion matrix M is similar to an upper triangular matrix. This 
can be shown just as in the complex case. Using elementary Gaussian operations, 
the general case can be reduced to the case of 2 x 2 matrices, where one wishes find 
a G HI so that 

★ *\ _ / 1 0\ fa b \ f 1 
-k) ~ \a l) \c d) \-a 1 

given a, 6, c, d G H. By permutation similarity we may assume that 6^0. Solvability 
for a is expressed by the non-commuting quadratic equation 

a 2 b + a(d — a) — c = , 

which always has a solution [^5[ ^6) . 

Note that in the complex case the similarity matrix obtained in this procedure is 
not in general unitary; however, a different procedure, Schur's lemma, triangularizes 
the matrix using unitary similarity. Schur's lemma has been extended to quaternionic 
matrices pp[ |. 

A modified version of the Jordan canonical form is valid for quaternionic matrices. 
Namely, every matrix M G H" x " is similar, over the quaternions, to a complex Jordan 



matrix J, defining a set of n complex eigenvalues. However the eigenvalues A^ 6 C 
are determined only up to complex conjugation [ fL8| . 

The Schur and Jordan canonical forms are associated with right eigenvalues 
Mif) = ipq, V 7 G HI " x 1 , gGi, which are determined only up to quaternionic similar- 
ity. This is further discussed in §, ||, §9|. 

Let us denote by Z[M] the complexification |Q, of the quaternionic 

matrix M , i.e. 

/ M —M* \ 

(1) Z[M] := I A £ M : J , M = M 1 + j M 2 , M M eC"". 

It has been shown in [jl8| that if J is the complex Jordan form of M then J © J* is the 
Jordan form of Z[M]. Consequently the spectrum of Z[M] is {A 1; A*, . . . , A„, A* }. 

4. On extending the determinant. The difficulty in extending the determi- 
nant to quaternions results from lack of commutativity. Starting with Cayley himself 
fl3f| all direct attempts at generalizing the concrete expression for the determinant 
have failed. Let us consider the case of 2 x 2 matrices. Here, one may consider four 
different generalizations: 

(2) ad — cb , ad — be , da — cb , da — be . 

It is easy to see that none of these expressions, alone or jointly, has relevance to the 
invertibility of the matrix. Consider, for example, the two matrices 

< 3 > A - A ( ) I ) ■ B - A ( )\ 

In the case of A, exactly two expressions in (||) vanish; in the case of B, all the four 
expressions are zero. However, both A and B are unitary. 

In a different line of attack, one may look for multiplicative functionals T from 
i* x " to H which coincide with the determinant on complex matrices. Again, the 
result is negative: 

There is no multiplicative functional 

T : M n * n — ► H , 

which coincides with det on complex diagonal matrices. 

It is enough to obtain one counterexample for n — 2. Consider the 2x2 matrices 

M=[ l V °A . N=( l V °. ), S {1 ° 



i J ' \ -i J ' V j 



Since S is invertible, J-[S] ^ 0, see Lemma |5.2| . Since SM = NS, we conclude 
that T[S] T[M] = T[N] T[S] , hence T[M] and F[N] should be similar. This is a 
contradiction because obviously Re{J r [M]} ^ Re{J r [A r ]}. 
□ 



and satisfies the scaling condition 

(6) V[ql] = \q\ n , Vq E 



-t— > 

<D 

X 

t-H 

5. On extending the absolute value of the determinant. A more positive 
result is obtained with respect to the functional |det|. 

Theorem 5.1. Sdet is the unique functional 

(4) V : H nxn 
which is multiplicative, i.e. 

(5) V[MN] = V[NM] = V[M] V[N] , 

Oh 

o 

Before proving this theorem, we make some observations concerning nonncgative 
multiplicative functionals. The only non-trivial part here is the last property, which 
has been proved elsewhere. 

LEMMA 5.2. If J- : H" x " — » M + is a non-constant multiplicative functional, 
then: 

1) F[S] = l, ifS 2 = I; 

O 

2) TiS^iS- 1 } = 1 and ^[S^MS] = T[M], if S is invertible; 
3] J r [i , l = 1 for all permutation matrix P; 

4) jF[M] = if and only if M is singular. 
Proof. 

> 

Multiplicativity and non-triviality implies that = 1. Now items 1-2 become 
trivial consequences of multiplicativity. Item 3 follows from the fact that every per- 
mutation matrix is a product of elementary permutation matrices Pi with Pf = I. As 
for item 4, if M is not singular then applying T to MM^ 1 — I implies that JF[Af] ^ 0. 
If M is singular, T ± by a result of §|, see @ pag. 58. 
□ 



Proof of Theorem 5. 1 



Let V be a functional satisfying (|]]|, [j)- Let {JSjj}™,- =1 be the usual canonical 
basis over HI in I" x ". Let g £ i be non-zero. Consider the n diagonal elementary 
matrices Mi(q) 



Mi{q) :=I+(q-l)E u 
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and the n(n — l)/2 upper triangular elementary matrices 



Mijiq) :=I + qE i0 , i <j . 

First we show that 

(7) F[Mi(q)] = \q\ • 

Indeed, by permutation similarity we see that !F[Mi(q)] is independent of 1 < i < n. 
So set f(q) := f[Mi(q)]. We have ql = n™ =1 Mi(q), hence 



T[qI]=\{T[Nh{q)]=P{<D 



i=l 

Hence f(q) — \q\, as claimed. Next we show that 

(8) nMij(q)] = 1 . 

Indeed, it is easy to see that 

Mr/(q) = M^-q) = M^-l^^M^-l), 

hence ^[M^ {q)\ = T[M l3 (q)}. Now (|) follows by multiplicativity. 

We have therefore established that T : H " x " — > H satisfies the three Dieudonne 
conditions (^), (||) and item 4 in Lemma [5.2| . Therefore, according to Dieudonne's 
result [jl6| = Ddct 21 = Sdct r for some r G M. Due to (||), it is readily seen that 
r = 1. 
□ 

Note that in general if T is multiplicative and r £ K then T r is also mul- 
tiplicative. Therefore, we have a one-parameter group of nonnegative multiplicative 
functionals: {Sdet r : r £ K.} (The case r = is interesting: it leads to the functional 
who se value is 1 on all the invertible matrices and otherwise). In view of Theorem 
(5.1) we conclude that this one-parameter family, plus the two constant functionals 
J-o[M] = and ^[M] = 1, are the only nonnegative multiplicative functionals on 
quaternionic matrices. 



6. Concrete description of the Study determinant. Theorem 5.1 has the 
following main corollaries: 

Corollary 6.1. If M is upper triangular then Sdet(M) = 11™= i \Mu \ . 

This follows easily by writing M explicitly as a product of elementary matrices, 
using (0|). 
□ 

Corollary 6.2. For all matrix M, Sdet(M) = n"=i M where Aj are the 
eigenvalues of M. 

By Lemma 5.2 item 2, it is enough to consider the Jordan form, or the Schur 



form, of M which is of the type considered by Corollary 3.1 
□ 



Since the eigenvalue identity just exhibited, restricted to complex matrices, is also 
valid for |det|, we get immediately: 

Corollary 6.3. For complex matrices we have Sdet(M) = |det(M)| . 

Let us define the adjoint of M by (M*)^ = M~. A matrix U G H™ x " is called 
unitary if U + U = I. According to the quaternionic Schur lemma p0| , every n x n 
quaternionic matrix, M, can be written as M = U + TU where U is unitary and T 
triangular. Since in addition Eq. (^) is obviously valid for lower as well as upper 
triangular matrices, we get: 

Corollary 6.4. Sdet(M+) = Sdet(M). In particular Sdet(f7) = 1 if U is 
unitary. 

The identity Sdet(M) = 1 may be taken as a basis to define the group of uni- 
modular matrices. 



Next, we calculate T in terms of singular values. The singular value decomposi- 
tion, SVD, for complex matrices extends to quaternionic matrices in a straightforward 
way. Every n x n quaternionic matrix, M, has the SVD M = UYjV where U and 
V are unitary, £ = Sj © 0, T, 1 = <j x © . . . © <Tk, where <j x > a 2 > . . . <Jk > are the 
singular values of M |3^, |3^, In these terms the following holds: 



Corollary 6.5. T[M] = JTJ 



7. Hermitian matrices. A quaternionic matrix H is called hermitian if H + = 
H . As we saw in section (Q), the common determinant cannot be extended to quater- 
nionic matrices. However, it can be extended to hermitian quaternionic matrices. The 
usual definition of the determinant in terms of permutations was generalized in the 
Chinese litera ture , see for example ref. |2^| . Another possible definition is analogous 
to Corollaries 3.2 and 6.5: 



\H\ 



Note that for hermitian matrices the eigenvalues are uniquely determined and real. 
This follows from the fact that Z[M] is also hermitian. Note that the set of hermitian 
matrices is not closed under products, and the functional det : H — » \H\ r is not 
multiplicative. However, it is invariant under congruence. 

It is easy to show that for hermitian matrices the following are equivalent: 
1) H is positive definite, i.e. x + Hx > for all non zero lei"* 1 ; 

2. All the eigenvalues Xi are positive; 

3. All the (signed) real determinants of the principal minors are positive. 

We conclude this section by comparing the functional Sdet[M], the functional 
\H\ r just defined, and the q-determinant fl7| 



\M\ q = det {Z[M]} , 



when Z [M] is defined in equation (ffl). From previous considerations, we have 

n 

\M\ q = ]J |Ai| 2 = Sdet[M + ] Sdet[M] = Sdet 2 [Af] = \M+M \ r . 

i=i 

Using this equation, one can extend the definition of Sdet from square to non square 
matrices. This approach is found in p^ ], where the resulting functional is called 
double determinant. 

8. Schur complements. Let 1Z be an associative ring. A matrix M £ TZ" x n 
is called invertible if MN — NM — I n for some N £ TZ " x " , which is necessarily 
unique. It is shown in (It]] that in case TZ = H, MN — I n implies NM — I n . 

The Schur complements procedure [ ^5| is a powerful tool in calculating inverses 
of matrices over rings. Let us write a generic ro-dimensional matrix M in block form 



M = 



A B 
C D 



Assuming that A £ 7Z k x k is invertible, one has 

h \ / A \ / I k A-'B 
CA- 1 I n - k Mo A s { I n - k 



(9) M = 

with 



A q := D - CA-'B 



We shall call A s the Schur complement of A in M. 

The invertibility of A ensures that the matrix M is invertible if and only if A s is 
invertible, and the inverse is given by 



(10) M- 1 = 



h -A- l B \f A- 1 \ / I k 
In-k Mo As- 1 \ -CA- 1 I n - k 



The inversion of an n-dimensional matrix is thus reduced to inversion of two smaller 
matrices, A £ lZ h x h and A s £ TZ { " _ fe) x ( ™ " fc) (plus some multiplications); repeated 
use of this size reduction can be used to invert the matrix efficiently. It is not as 
efficient as Gaussian elimination, but the latter may not be available in general rings. 



Corollary 8.1. Sdet ^ ^ ^ J = Sdet [A] Sdet[D - CA- X B\ as long as 
A" 1 exists. 

Indeed, from the construction of Sdet in the last section, we see that its value on 
each of the two block-triangular matrices in (||) is 1; since the eigenvalues of a direct 
sum are the union of the eigenvalues of the summands, we get that Sdet [A © As] = 
Sdet [A] Sdet [As]. This plus multiplicativity implies the result. 
□ 







[(c 


2)] 



As a result of the Schur complements determinant formula just exhibited, we get 
the following commutation formula for Sdet, which generalizes a well known property 
of det (actually, of |det|): 



X 

in 



Corollary 8.2. Sdet [J + MN] = Sdet [7 + ATM] for all M G 
7VeH mx ". 

Indeed, consider the matrix 



and 



o 
o 

o 
o 

> 

O 

(D 

O 

OH 

o 
-a 

> 

O 

<D 

o 



M I, 



and apply to it Schur complements with respect to both I x and I 2l respectively. We 
get 



and 



Sdet 



Sdet 



\( h 




\ M 


D\ 


\( h 




\ M 





Sdet [/J Sdet[/ 2 - MI[ X N] , 



Sdet[7 2 ] Sdet[/ 2 - NI^M] , 



implying the identity. 

□ 



9. The case of 2 x 2 matrices. In this last section, we discuss inversion, adjoint 
and determinant for 2x2 quaternionic matrices. 

• Inversion. Let 



M = 



a b 
c d 



be an invertible 2x2 matrix with quaternionic entries. When a, 6, c, d are all non-zero, 
four parallel applications of the Schur complement formula (|lo| ) lead to a concrete 
description of the inverse: 



(11) 



where 



(12) 



M" 1 



a b 
c d 



d = (a — bd 1 c) 1 , 
5= (b-ac-'d)- 1 



b=(c-db- 1 a)- 1 , 
d= {d-ca-'b)- 1 , 



see Giirsey j3<| page 115. The invertibility of M guarantees that these four values 
are well-defined non-zero quaternions. What happens if some of the entries of M 
vanish? Assume for example that a = 0. The invertibility of M implies that b, c ^ 0. 
Consequently, the element d — ca~ 1 b has infinite modulus. In this case, we define 



d := lim (d — ca 1 b) 1 . 

a—>0 



A simple calculation, 



\d\ := lim r— -r = lira -r— n ; -r 

a^o Id-ca- 1 ^! q^o \c\ \c- 1 d — a -1 ©! 



= um m ; — rr = um , , ,, , = , 

a^o |c| \ac- x d-b\ a^O |c| |6| 

shows that d = 0. Thus, 



We conclude that Eq. (|12|) remain valid under appropriate conventions, when some 
entries in M are zero. We do not have a clear generalization of this phenomenon for 
n > 2. 

• Adjoint. Eq s. |l2] ) are valid in every associative ring 1Z. In case 1Z is also 
commutative, Eqs. ( |1 1|J12| ) reduce to the well known formula 

(13) M" 1 = . 
1 ; det[M] 

In calculating the inverse of real and complex matrices, (|l^) is of great theoretical 
importance. So far, we have failed to generalize this formula to quaternion matrices. 
At first sight, it might make sense to conjecture a non-commuting expression of the 
general form, 

(14) M" 1 = PAdj[M] Q , 

with quaternionic diagonal matrices P = diag{p!,p 2 } and Q — di&g {q 1 , q 2 } . Never- 
theless, the resulting constraints 

Pi = aq^- 1 , p 2 = -cq^- 1 , 

Pi = ~bq 2 b- 1 , p 2 = dq^- 1 , 

which, for commutative fields, are satisfied if P = det^ x [Af] / and Q = I, are not al- 
ways solvable. For example, the first matrix in cannot be written in the form (fL4|). 
Whether a further weakening, beyond (|l4|), of formula ( |l3| ) is valid for quaternionic 
matrices remains an open problem. The mere definition of Adj[M], M 6 H" x ", 
n > 2, preserving (|l3|), is not clear. 

A different generalization of ( |l3| ) for 2x2 quaternionic matrices may be obtained 
using a Hadamard product between a non negative matrix and a termwise-unitary 
quaternionic matrix 



(15) M 



-i 1 ( \d\ \b\ 



Sdet[M] V l c l 1° 




Another description of the inverse matrix is offered in equation (37) of the paper 
of Chen il. 
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• Determinant. For n = 2, it is noteworthy that the following four quaternion 
expressions are equal: 

\a\ \d - ca- x b\ = \b\ \c - db^a] = |o| \b - ac^d] = \d\ \a - 6d _1 c| . 



.1 



it follows that each of this 



From the Schur complement formula, Corollary 
expressions, properly extended in case a, b, c or d is zero, expresses the value of 

Sdet''" 



Applying this formula on the two unitary matrices in (^|), one obtains 

the expected result (these matrices are unitary, hence unimodular). For hermitian 
quaternionic matrices, the real determinant is given by 



/ a 




V 9 


s) 



Ai A 2 = aS - \q\ 



a , S e 



q e 
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